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Abstract: We have developed Water Vapour Radiometers (WVRs) for the Australia Telescope Com- 
pact Array (ATCA) that are capable of determining path fluctuations by virtue of measuring small 
temperature fluctuations in the atmosphere using the 22.2 GHz water vapour line for each of the six 
antennae. By measuring the line of sight variations of the water vapour, the induced path excess and 
thus the phase delay can be estimated and corrections can then be applied during data reduction. This 
reduces decorrelation of the source signal. We demonstrate how this recovers the telescope's efficiency 
and image quality as well as how this improves the telescope's ability to use longer baselines at higher 
frequencies, thereby resulting in higher spatial resolution. A description of the WVR hardware design, 
their calibration and water vapour retrieval mechanism is given. 
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1 Introduction 

1.1 Overview 

The Australia Telescope Compact Array (ATCA) is a 
radio interferometer capable of observing in the cen- 
timetre and millimetre wavelength regime of the elec- 
tromagnetic spectrum. It was for a little over a decade 
the only millimetre interferometer in the southern hemi- 
sphere. 

Variations in the water vapour content in the atmo- 
sphere along the lines of sight from the source being 
observed to each antenna lead to fluctuations in the re- 
fractive index, and so to path delays and phase fluctu- 
ations. At millimetre wavelengths, these fluctuations 
are rapid, requiring frequent measurement of a cali- 
brator source in order to be followed. This may be as 
often as every 5-10 minutes, depending on weather con- 
ditions and frequency. This can consume up to a quar- 
ter of the observing time. The ATCA has 3 mm re- 
ceivers (tuneable from 83 to 105 GHz), 7 mm receivers 
(tuneable from 30 to 50 GHz) and 15 mm receivers 
(tuneable from 16 to 26 GHz). They are mounted 
on a turret which moves the feed for the requested 
receiver on axis automatically. The water vapour in- 
duced phase variations increase with the length of the 
baseline between antennae. With the Australia Tele- 
scope Compact Array, the largest useable baseline at 3 
mm is found to be ~ 300m, yet the maximum baseline 
separation of the five antennae equipped with 3 mm 
receivers is 3km (twice that, 6 km, for 7 mm and 15 
mm), and this is used regularly at centimetre wave- 
bands. If the maximum available baselines could be 
used at millimetre wavelengths it would result in up to 
a factor of ten improvement in the spatial resolution 
for millimetre imaging at 3 mm. On shorter baselines 



improved phase stability would also improve the signal 
to noise ratio of millimetre-wave observations. 

It may also extend the millimetre wave observing, al- 
lowing the telescope to be used in a wider range of 
conditions. In 2006, 33% of proposals in the April 
semester involved observations at 3mm: in 2007 the 
figure was 25%. In 2008, 37% of Australia Telescope 
Compact Array proposals requested 3mm observing, 
and fully 45% requested 3mm and/or 7mm observa- 
tions. These figures are comparable throughout the 
years up until the April 2011 semester. The interest in 
millimetre observing has remained strong with the in- 
stallation and commissioning of the Co mpact Array 
Broad band Backend (CABB) in 2009 ([Wilson et all 

(HqU)). 



1.2 The Atmosphere 

The various constituents of the atmosphere all have 
an effect on the transmission of signals from space. In 
the millimetre region (between about 16 GHz to 300 
GHz), the main contributors to atmospheric attenu- 
ation are O2 and H2O. For A < 1.5 cm, the optical 
depth of the troposphere becomes significant, lead- 
ing to increased system temp e rature s due to atmo- 
spheric emission l|Carilli et all |l999)). The refrac- 
tivity of water vapour at these frequencies is about 
20 times greater than in the ne ar-infrared or optical 
regimes (|Thompson et all l|200lM . Figure [T] shows the 
sky brightness caused by these opacity contributions. 
The water vapour lines at 22.2 and 183.3 GHz are 
prominent. 
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When evaluating the atmospheric conditions, several 
models were used to characterise the temperature dif- 
ferences arising due to variations in the water vapour 
content along the lines of sight. The first was the 
Millimeter P ropagation Model (MPM) developed by 
iLiebd l|l985T ) and implemented in the Miriad software 
package. While accurate enough for purposes of fil- 
ter op_tiinisa^ionj_the model has been improved upon 
bv lPardo et al.l lj200lT ) who developed the Atmospheric 
Transmission at Microwaves (ATM) model used to sim- 
ulate PWV induced excess temperatures. The well 
characterised VanVleck - Weisskopf (WW) profiles 
were also used to determine the model atmosphere pa- 
rameters for the water vapour calibration factors as 
well as for the weighting coefficients determined in H4.2I 
to combine measurements made at several frequencies 
across the 22.2 GHz water line. The same WW pro- 
files used for the coefficients were also incorporated 
into the ATM model. 



Atmospheric Emission Spectrum frem 1 GHz te 200 GH; for Ele*atlcn = 2c5m P = lQ00hPo T=29DK 
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Figure 1: The atmospheric emission from 1 GHz 
to 200 GHz for a range of PWV co nditions, mod- 
elled using the iPardo et al.l (|200ll ) ATM model. 
Variations in precipitable water vapour (PWV) are 
shown from 1 mm to 30 mm with a constant atmo- 
sphere at 235 m elevation, 1000 hPa pressure and 
a temperature of 290 K. The features at 22.2 GHz 
and 183.3 GHz are the prominent water vapour 
lines. The 60 GHz and 118 GHz features are 
caused by O2. It is immediately evident that for 
the PWV values of 5 mm to 30 mm encountered at 
the ATCA site in Narrabri, only the 22.3 GHz line 
is a viable candidate as the 183.3 GHz line is com- 
pletely saturated. For locations where the PWV 
falls below 2 mm, the 183.3 GHz line is the better 
choice. Note also that the continuum increases as 
r; j/ 2 , arising from the liquid contribution. 



2 WVR Design 
2.1 Design principle 

The WVRs developed for ATCA are an uncooled de- 
sign based on four filters placed between the frequen- 
cies of 16 - 26 GHz, each with 1 GHz of bandwidth 
to increase the sensitivity. This allows the WVRs to 



measure small variations in the spectrum of the 22.2 
GHz water vapour emission and derive phase vari- 
ations based on these measurements. They are co- 
located with the science receivers, with the WVR feeds 
offset by 13 cm to the on axis mm receiver being used, 
hence look through almost the same beam onto the 
sky. They operate at frequencies the antennae were 
designed to perform at. They further are low mainte- 
nance due to the lack of a cryogenically cooled stage. 
By using a differential water vapour determination, ab- 
solute measurements are not required. We use them 
to measure the minute temperature differences caused 
by the fluctuations in the water vapour between the 
lines of sight for each antenna pair. 

2.2 Determining a figure of merit 

We introduce a figure of merit to model performance 
of the WVRs in order to determine the accuracy to 
which the temperature fluctuations must be measured. 
The expected level of decorrelation of the incoming sig- 
nal due to atmosp heric effects can be expressed using 
the Ruze formula (jRuzd (|1966T shown in equation [T] 
This equation states that for a signal reflected off an 
antenna surface of area Aq, the aperture efficiency is 
A/A . 

A = A e~^ 2 (1) 

where a is the RMS surface accuracy of the an- 
tenna. Applying this to interferometry, the phase shift 
inducing error term which for single dish application 
is Ana/X is halved to 2ira/X because the locus of the 
error is above the antenna and the path error is thus 
induced only once, a is now the path difference caused 
by the water vapour fluctuations. By expressing the 
required precision N as a fraction of the wavelength 
A, i.e. 1/N — u/X, the figure of merit is obtained, 
hereafter referred to as the correlation efficiency e: 

e = e -(^) 2 ( 2 ) 

Equation [2] is plotted in Figure [2] to show the A 
fraction to which one needs to correct in order for the 
WVRs to provide useful levels of improvement: For 
a correlation efficiency of at least e = 0.5, the water 
vapour induced path delay must be determined to bet- 
ter than Nthreshoid = A/7.5. 

2.3 Filter optimisation 

The original filter locations were based on a design 
by Robert J. Sault. The frequencies were chosen to 
cover the water line in its most sensitive regions based 
on atmospheric modelling and a 10 year archive of ra- 
dio sonde data providing the input for the model. We 
have repeated these optimisation steps but modified 
them for a wider filter bandwidth with 3 filters placed 
across the line and one filter on the continuum. We 
have also identified a number of RFI sources (chiefly 
communication satellites) and have managed to keep 
the filter locations outside the regions of the commer- 
cially allocated radio spectrum. The filter placement 
is illustrated in Figure [3] Shown are the overall wa- 
ter line excess variations for a 10 year period using 
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Water Vapour Line Excess Temperatur 



A fraction (N) vs. Correlation efficiency (e) 
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Figure 2: Fraction of wavelength vs. the result- 
ing correlation efficiency: The solid line plots the 
Ruze formula (see equation Q]) . Over-plotted are 
shown the resulting correlation efficiency from the 
measured phase noise in test observations. Stars 
show short baselines with interpolated data, where 
correlation efficiency is generally better than 0.95. 
Plus signs show these baselines with WVR correc- 
tions applied. While this results in added noise, 
the lowered correlation efficiency remains above 
0.9. On the long baselines, the WVR improve- 
ments are substantial: without corrections, effi- 
ciencies are about 0.65 (squares), with corrections 
better than 0.9 (diamonds). 



radio sonde data as input for the model. The dotted 
lines represent to maxima/minima values encountered 
in the entire data set. 



2.4 Determining Noise Limits 

The theoretical sensitivity of an ideal total-power ra- 
diometer with no gain fluctuations is calculated as 
(|Ulabv et all (|l98lh h 

where AT is the minimum detectable change in the 
radiometric antenna temperature T a , T e is the noise 
temperature of the entire signal chain, r is the inte- 
gration time and A/ is the bandwidth of the system. 
This is illustrated in Figures [3] and [5] Solutions both 
for the sensitivity Q as well as the minimum time re- 
quired to reach the target sensitivity of 12 mK ((5JI are 
given. 

T e , the total noise temperature, is calculated as a 
cascaded system taking each component's noise tem- 
perature contribution in the s ignal c h ain in to account 
as given by the Friis formula (|Krauj i|l966D ): 



0.25 — 




Figure 3: Filter placement over the 22.2 GHz wa- 
ter line: Shown are the overall water line excess 
variations for a 10 year period using radio sonde 
data as input for the model. The dotted lines rep- 
resent to maxima/minima values encountered in 
the entire data set. 



T ei + 
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E 

i 



Te, 



(4) 



where T ei is the noise temperature of the first stage, 
T ejv and Gn are the noise temperature and the gain 
of the N th stage respectively. 

To specify the noise limitation, the fundamental equa- 
tion determining a radio receiver noise floor is appli- 
cable. The power received is determined by: 



P = fcTA/ 



(5) 



where k is the Boltzmann constant, T is the tem- 
perature in K and A/ is the bandwidth of the receiver 
in Hz. Note that in high frequency engineering, it is 
useful to express the power and noise ratios as refer- 
enced to 1 milliwatt (abbreviated as dBm). P hence 
is converted to mW for all calculations. Assuming 290 
K for the temperature and 1 GHz of bandwidth, the 
minimum equivalent input noise of the water vapour 
radiometer is -83 dBm, also noted on the receiver di- 
agram in Figure [6] The noise floor Nj can then be 
determined: 



N f = -174dBm + iV /i9 + 10 



(A/) [dBm] (6) 



where Nf ig is the noise figure of the system and 
-174 dBm corresponds to the input noise of 1 Hz of 
bandwidth at room temperature (290 K). The noise 
figure, noise floor and receiver temperatures have been 
determined for all channels and units and are listed in 
Table rj] 

2.5 Hardware 

The uncooled radiometers are designed to reside inside 
the vertex room, attached to the millimetre receiver 
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Antenna 


Unit 


Filter [GHz] 


Tree [K] 


Noise Figure [dBl 

O LI 

Nfig 


Noise Floor [dBm] 
N f 


1 


1 


16.5 
18.9 
22.9 
25.5 
TP 


319.2 
251.4 
249.9 
316.0 
422.4 


3.22 
2.71 
2.70 
3.20 
3.20 


-81.43 
-81.84 
-81.30 
-80.76 
-70.80 


2 


2 


16.5 
18.9 
22.9 
25.5 
TP 


370.3 
266.8 
266.8 
334.6 
438.3 


3.57 
2.83 
2.83 
3.33 
3.33 


-81.13 
-81.48 
-81.17 
-80.46 
-70.67 


3 


3 


16.5 
18.9 

22.9 

op; k 
zo.d 

TP 


300.8 
267.0 
287.0 

386.1 


3.09 
2.83 
2.99 

o.oU 

3.30 


-81.51 
-81.48 
-80.76 

— oU.o ( 

-70.70 


4 


7 


16.5 
18.9 
22.9 
25.5 
TP 


305.0 
251.8 
320.5 
382.5 
398.1 


3.12 
2.71 
3.23 
3.65 
3.65 


-81.58 
-81.70 
-80.72 
-80.09 
-70.35 


5 


5 


16.5 
18.9 
22.9 
25.5 
TP 


288.8 
279.3 
316.8 
352.5 
620.8 


3.00 
2.93 
3.21 
3.45 
3.45 


-81.81 
-81.48 
-80.67 
-80.55 
-70.55 


6 


6 


16.5 
18.9 
22.9 
25.5 
TP 


296.8 
244.2 
265.6 
295.0 
676.5 


3.06 
2.65 
2.82 
3.05 
3.05 


-81.64 
-81.76 
-81.01 
-80.95 
-70.95 



Table 1: The receiver temperatures T rec , the noise figure and noise floor for each filter in each of the 
WVRs. TP is the total power channel, which is 10 GHz wide, from 16 to 26 GHz. Note that Unit 4 (for 
historical reasons) remained with the manufacturer for extended testing purposes and antenna 4 thus was 
provided with unit 7 instead. 
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Integration Time vs. Sensitivity 



T a +T e vs. Time to reach 12mK Sensitivity 
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Figure 4: The theoretical sensitivity of a radiome- 
ter with 1 GHz bandwidth as given by equation [31 
calculated for summed antenna and noise temper- 
atures of 290 K, 400 K, 600 K, 800 K and 1600 K 
respectively. 



system mounting cage, also referred to as the "mil- 
limetre package". The WVR feed horn is situated in 
such a way that regardless of which mm receiver is on 
axis, the WVR feed is offse t by no more t han 1 2 cm 
from the axis. According to lCooper et al. (1992), the 
ATCA antennae have a diameter of 22 m, and an / 
ratio of 2.0, so the beams for the receiver and WVR 
are separated only by 5.4 m at 2000 m altitude, thus 
still overlapping by 68% of the beam footprint. The 
majority of troposphere induced path delays originate 
from below that altitude. Hence the WVRs are mea- 
suring fluctuations through virtually the same column 
of atmosphere as the science receivers. The beams are 
completely separated above an elevation of ~ 8,000 m. 

To illustrate the signal path and the components in- 
volved in the water vapour radiometers, as well as to 
visualise the insulation characteristics of the package, 
we now walk through the signal path while simultane- 
ously looking at the following two figures: 

• The radiometer schematics in Figure [6l starting 
at the top, 

• Unit 4 in Figure [Jj 

The signal path is entered at the feed horn (only 
shown on the schematic) where a weak signal of -83 
dBm at 1 GHz of bandwidth (corresponding to the 
widest filters in the system) is received. The waveg- 
uide (numbered 1) feeds this signal into the first two 
low noise Miteq amplifiers (LNA's, numbered 2) which 
amplify the signal each by +36 dB. The signal contin- 
ues via an Inmet matching pad through a 26 GHz low 
pass filter (numbered 3) and another matching pad, in- 
troducing a loss of -25 dB. Then the first MCLI splitter 



Figure 5: The integration time required to reach 
a sensitivity of 12 mK as a function of the sum of 
antenna and noise temperatures T a + T e terms. 



is reached (numbered 4) which splits off the total power 
signal from the other filters. The total power channel 
now directly feeds into the Herotek tunnel diode de- 
tector and from there into the Analog Devices ultra 
precision operational amplifiers from where the signal 
is digitised and fed into the dataset (this last step is 
not shown). Following the other signal path after the 
total power splitter, another matching pad is entered 
and another Miteq LNA (numbered 5) amplifies the 
signal by +36 dB, then the signal is fed through a 
Ditom isolator, an MCLI four-way splitter (numbered 
6) and into the ASTROWAVE band pass filters (num- 
bered 7), one for each of the filter frequencies of 16.5, 
18.9, 22.9 and 25.5 GHz respectively. Then the sig- 
nal for each filter follows the Herotek tunnel diode de- 
tectors into the amplifiers and then into the datasets 
where the data is obtained via the mm receiver pack- 
age data bus. The datasets are providing digitisation 
of the analog signal voltages from the tunnel diodes. 
The arrangement of the WVR feed as well as the hard- 
ware locations are shown in Figure [9] showing a view 
from the top of the focal plane into the feed horns sit- 
uated on top of the mm package dewar, as well as in 
Figure [10] the WVR unit can be seen mounted on the 
side of the mm receiver package. 

The required temperature stability was derived from 
the overall ga in stability requirements discussed by 
iBremerl (|l997t ). The most temperature dependent com- 
ponents in the WVR are the RF amplifiers which have 
a gain change of about 0.04 dB per 1 K. This implies 
a temperature control requirement of ± 1 mK at the 
RF plate control point. The limits were verified by 
calculating the Allan variance and deviation shown in 
equations [JJ and [5] This is a useful metric because sim- 
ilarly to a power spectrum it highlights behaviour in 
the time domain such as externally caused tempera- 
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Figure 6: A schematic diagram of the WVR kindly 
provided by Christoph Brem (CASS Narrabri). 
The dB numbers on the lcfthand side show the 
signal loss (or gain) at each stage along the signal 
path. The blue figures to the right are used as 
stage references in the text. 




Figure 7: A photograph of the RF plate of Unit 
4. The foam insulating material is visible around 
the unit and the components can be identified by 
comparing to Figure [6] 

ture oscillations e.g. by the vertex room air condition- 
ing system. The Allan deviation is calculated as the 
square root of the Allan variance, which in turn is one 
half the time average of the squares of the differences 
between successive data points of the frequency devi- 
ation, sampled over the sampling period as shown in 
equation 

4M = i<(TVTT-TV) 2 ) (7) 
and from this the Allan deviation: 

*t(t) = sftUr) (8) 

The RF Plate temperature control point data was 
found to be well within the required limits at values 
between 0.1 and 0.4 mK. One example plot is shown 
in Figure fTTT 

To minimise temperature changes with ambient 
temperature variation, multi layered insulation sur- 
rounding the RF plate was installed. The ideal shell 
material is anisotropic with good insulation 'across' 
the medium and good conductivity 'along' the medium 
(i.e. around the RF plate). A useful approximation to 
this is several alternate layers of thermal insulator and 
conductor. This presents a challenge if easy access to 
the RF plate is to be maintained. A practical solution 
was to build two insulating layers and three conductive 
(aluminium) layers with the inner two actively temper- 
ature controlled. The minimum gradient change was 
successfully achieved using this method. Figure [5] il- 
lustrates in a cross section how this was achieved: The 
top halves of the conducting layers were rendered as 
a transparent plastic in this visualisation to illustrate 
the layering. The hollow spaces in between are foam 
filled in the real units. 

For both gain stability and temperature stability the 
+ 10V drive voltage needs to be exceedingly stable, be- 
cause changes in power dissipation of RF components 
would change local temperature. Great care was taken 
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Figure 8: Cross section render through one of the 
WVRs. In the centre is the RF plate with com- 
ponents visible. The RF plate itself is shielded in- 
side a first aluminium enclosure, which itself is sur- 
rounded by foam (not shown in the render). Then 
a second aluminium shell surrounds this with more 
foam on the outside and finally the outside shell 
which is not thermally controlled. The outside is 
however painted white to minimise thermal energy 
uptake through solar radiation when pointing the 
antennas near the Sun . Visualised by the author 
using Maxwell Render ( Nextlimitl (|201ll )) based on 
an AutoCAD model. 



to stabilise this using precision voltage references in- 
stalled on the RF plate PCB. This helped maintain 
the temperature stabilised environment with all of the 
control semi-conductors (which all have varying power 
dissipation) outside of the thermal enclosure directly 
on the main heat sink. 

To thermally isolate the RF plate from the environ- 
ment, all electrical conductors leading into the RF elec- 
tronics were fabricated as thin traces on a 200mm long 
specially manufactured (and flexible) Kapton printed 
circuit. The input coaxial feed is fabricated from stain- 
less steel for low thermal conductivity, with thin silver 
layers on the centre conductor and inside of the outer 
conductor for low electrical loss. This reduces thermal 
leakage but it is still the weakest link in the thermal 
design because this cable penetrates the shield. 

Sharp cut off filters to minimise overlap between the 
two closest channels were needed. A 7th order Cheby- 
chev filter design was chosen as they give the sharpest 
cut off and the theoretical f dB ripple in the pass band 
is dampened out by the insertion loss. 7th order also 
gives a rather high insertion loss (5 to 7 dB) but this 
is at high signal level at the end of the signal channel 
and there is sufficient gain available to render the loss 
minor. 

The filters were implemented as Micro-Strip filters etched 
on 0.25 inch backed 'Duroid' (i.e. ceramic loaded PTFEQ) 




Figure 9: Top view of the millimetre dewar with all 
the feeds visible. Clockwise from the top: WVR 
feed with gold coloured waveguide leading the sky 
signal to the WVR box mounted on the right hand 
side of the mm package. Next is the 7 mm feed, 
then follows 15 mm and lastly to the left is the 
3mm feed. 



substrate. This gave a small physical size, helping 
keep the overall WVR package small. Also, by keep- 
ing the RF plate as small as possible, the temperature 
gradients induced by a given thermal leakage remain 
smaller. Further, this makes for a mechanically very 
solid, stable, filter. 

There was some concern regarding thermal stability of 
this substrate material but bench tests of a completed 
filter over a 50 K range indicated that the insertion 
loss had a temperature coefficient of only 0.01 dB / K, 
or 0.23 parts in 10 4 / 10 mK, which is well below the 
RF amplifier's temperature coefficient of 0.04 dB / K. 



2.6 Calibration 

Calibration of the water vapour radiometers was achieved 
by measuring the individual channel voltage at two 
known set points. The first set point measures the hot 
load (or paddle), consisting of activating the paddle 
mechanism and reading the temperature of a probe 
inside the absorber foam that is mounted on the pad- 
dle. There is also a case temperature probe (called T5) 
on each water vapour radiometer which can be used to 
indicate the ambient temperature (and thus approxi- 
mates the hot load temperature). The colder set point 
is achieved by setting a styrofoam box on top of the 
feed horn with a microwave absorber inside submerged 
in a sufficient amount of liquid nitrogen so that the ab- 
sorber along with the LN2 becomes optically thick. 

With cold and hot load on the feed, the spillover tem- 
perature is not a factor and the gain G can be deter- 
mined via the Y factor: 



ipTFE = Polytetrafluoroethylene, more commonly 
known as the brand name Teflon. 



Tuot + Tr c 

Tcold + TR e 



VHot 
Vcola 



: Y 



(9) 
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Figure 10: A side view of the millimetre package 
with the WVR mounted on the side (the white 
box). The WVR paddle, feed horn as well as 
waveguide can be clearly seen. The WVR pad- 
dle is in the same position as in Figure [9l i.e. not 
obstructing any of the feed horns. Photo by Peter 
Mirtschin in April 2011. 

where Tuot is the hot load temperature, Tcoid is 
the cold load temperature and 

„ Tnot — YTcoid , 1fV . 

iRec = - (10) 

Vffot and Vcoid are the hot and cold load voltages 
measured and 

f~, _ Tuot + Tfl ec . . 

_ v \ ' 

VHot 

The gains and Y factors have been determined 
three times, the data for unit 7 (antenna 4) and its 
first two calibrations are shown in Table [2] The Y fac- 
tor is seen to change by less than 1% between the cal- 
ibrations which were spaced 3 months apart. A third 
calibration was done 7 months after the first calibra- 
tion and found no significant change in Y factors. This 
means that with the hot load alone, gain drift correc- 
ti ons can be applied. This is discussed in further detail 
in llndermuehlell|201ll ). 

After each calibration series, a zenith sky spectrum 
was obtained for each antenna. Figure [T2] shows this 
spectrum for Unit 1. The sky temperatures in the 
channels follow the model atmosphere as expected and 



Allan Deviation of RF Plate Temperature for Unit 1 











\ - 















100 1000 toooo 



Figure 11: The Allan deviation of the RF Plate 
temperature control point for Unit 1 in the evening 
of August 23 2011. The dotted lines represent 
the maximum temperature drift as per the de- 
sign specifications. The temperatures therefore are 
more stable than the design criteria required by 
about 3 orders of magnitude for unit 1. 

outline the shape of the H2O line on top of the v 2 
continuum. The sky temperature in the total power 
channel (spanning the full 16 - 26 GHz) is at a simi- 
lar level to the highest continuum point at 25.5 GHz, 
and the shape of the water line is traced in all filter 
positions. A further representation of this data for all 
units is shown in Figure [13] which also shows the ATM 
derived model atmospheres for 5 mm, 50 mm and 23 
mm PWV reflecting the site minimum, site maximum 
and current conditions respectively when these mea- 
surements were taken. 

Unit 1 Zenith Spectrum 
50 r — 1 1 1 1 1 1 1 — 1 1 1 1 1 1 1 1 1 — 1 1 1 — : 

: : 

40 T O T Fll , er /' \ - 

: A T Totol Power /' \ ■ 




10 - 0--'"" — 

ob , , , I , , , I , , , 1 , , , 1 , , , : 

16 IE 20 22 24 26 

F,„„,nt, [GHz] 

Figure 12: Zenith sky spectrum for unit 1. The 
temperature excess reaching above the v 2 contin- 
uum of the three outside channels is approximately 
20 K. Data taken on uMarch 7 2011 in clear skies. 



2.7 Sky dips 

During the calibration runs, several sky dips were ob- 
tained, one set of which is depicted in Figure [141 These 
are measurements of the sky flux obtained at a num- 
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16.5 


18.9 


22.9 


25.5 


Total 




GHz 


GHz 


GHz 


GHz 


Power 


Y Factor Aug 2010 


1.582 


1.682 


1.571 


1.506 


1.490 


Y Factor Nov 2010 


1.579 


1.677 


1.570 


1.500 


1.476 


Y Factor difference 


+0.2% 


+0.3% 


+0.06% 


+0.4% 


+0.9% 


T Rec K Aug 2010 


304.1 


247.9 


311.4 


360.9 


375.4 


T Rec K Nov 2010 


306.0 


250.3 


311.9 


366.6 


388.3 


T Rec difference 


1.0% 


1.0% 


1.0% 


-1.5% 


-3.4% 



Table 2: Gain comparisons for calibrations executed on Unit 7 in August and November 2010. 



ber of elevation angles from 90 degrees to 13 degrees 
in steps of 10 degrees. Based on the skydip data, each 
channel for each WVR unit had its r, zenith Ts ys as 
well as its spillover term Ts determined by a least 
squares fit of equation [12] to the data. 

T secm = T S + T A (1 - e C— »)) (12) 
which, for small r becomes: 

T aec{6) =T s + T A {Tsec6) (13) 
where Ta is the atmospheric temperature. 

The temperature intercept derived in eg uation l 1 21 yields 
the spillover temperature. In the skydip data pre- 
sented in Figure Q3] we show the measured points at 
each elevation, connected with a dashed line. The in- 
tersect of the fit with the X axis yields the spillover 
temperature. The solid lines depict a linear fit through 
the skydip points. The fitted data to the skydip equa- 
tion [T^] is shown by dotted lines. Figure [Tol shows the 
skydip functions for all units and illustrates how the 
spillover temperatures are derived by intersecting the 
fit with the X axis. Table provides the spillover tem- 
peratures. We find that they are negative, i.e. the 
WVRs either underestimate the measured sky temper- 
atures near zenith, or they overestimate the sky tem- 
peratures near the horizon. This is greatest for the 
22.9 GHz filter in the water vapour line itself, and less 
for the other filters. These are, however, a small frac- 
tion of the ground temperature of ~290 K and amount 
to between 1% and 7%. It is also possible that the 
microwave absorber submerged in liquid nitrogen was 
warmer than the theoretical 77 K, either due to water 
freezing to the outside of the styrofo am box, LN 2 b e- 
ing warmer than the 77 K assumed (|Ericksonl (|l985l )) 
or a combination thereof. This then would result in 
negative values derived for Ts; the true values would 
be obtained by adding the error in T co id (if known). 

The variation in receiver temperature, which is as- 
sumed constant in this analysis, is also of order 1 - 
3.5% (see Table [2]). However, this result also indicates 
limitations in the simple sec((9) plus constant atmo- 
sphere model assumed to apply to the skydips. For a 
real atmosphere, the temperature will vary as a func- 
tion of height. The opacity is determined from the 
curvature in each skydip. The results from the fit, in- 
cluding the errors, are listed in Table [3] Further exam- 
ining the optical depths in Figure \W\ we see that the 



derived zenith opacities are well behaved when com- 
pared to the zenith opacities of the model atmosphere. 
The data is consistent with typical atmospheric con- 
ditions but the variations in the values determined 
also indicate that WVR data cannot be used for abso- 
lute determination of path lengths through the atmo- 
sphere. Only the differential path lengths obtained by 
subtracting the offsets between each antenna can be 
derived. 

2.8 Time domain analysis of the fil- 
ter signal 

A wavelet analysis has also been applied to the input 
signals received by the WVR units in order to verify 
that the signals they record are indeed dominated by 
the emission from the sky, and that they are tracking 
variations in the sky emission, as opposed to locally- 
based signals in the antenna. Wavelet analysis pro- 
vides a mean of examining the variation of power re- 
ceive d within a time series (e.g. iTorrence and Compol 
1998). This allows determination of not only the fre- 
quencies on which the dominant power is received by 
a system (as, e.g. a Fourier analysis would show), but 
also the time scales on which that power itself varies. 
We have used a Morlet transform for this analysis 
(which is equivalent to a plane wave modulated by 
a Gaussian) using IDL-based routine^] to process the 
input time series. We show here an analysis based on 
the signal received through the 25.5 GHz filter on each 
WVR unit, as well as, for comparison, the signal from 
the temperature sensor on the outer shell of each WVR 
unit. The latter is tracking variations in the tempera- 
ture of the receiver cabin (or vertex room). 

Figure [T7] shows the wavelet analysis for the raw 
voltages measured through the 25.5 GHz filter for a 
12 hr period on 23 May 2012 for each of the 6 anten- 
nae. The other filters show a similar response. Time 
is plotted on the z-axis and period on the y-axis and 
the intensity of the Morlet transform indicates where 
the power in the input signal arises in this domain. 
To first order all WVR units show a similar pattern. 
A series of weather systems passed across the obser- 
vatory over the second half of the day, for instance 
a weather front at ~15:45UT, and these are clearly 
tracked by each WVR unit, with substantial power 
fluctuations evident. These short-lived events gener- 

2 See http://paos.colorado.edu/research/wavelets/ 
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Antenna Filter [GHz] T S [K] T s ,< 





16.5 


-2.4 


0.50 


0.04 


0.001 


CA01 


18.9 


-3.7 


0.57 


0.07 


0.001 


22.9 


-6.5 


0.95 


0.20 


0.003 




25.5 


-5.6 


0.83 


0.13 


0.002 





16.5 


-6.9 


0.48 


0.04 


0.001 


CA02 


18.9 


-6.6 


0.56 


0.07 


0.002 


22.9 


-12.1 


1.60 


0.21 


0.005 




25.5 


-11.3 


1.14 


0.13 


0.004 





16.5 


-7.0 


0.47 


0.05 


0.001 


CA03 


18.9 


-9.0 


0.84 


0.08 


0.002 


22.9 


-14.1 


2.26 


0.23 


0.008 




25.5 


-8.8 


1.37 


0.13 


0.004 





16.5 


-10.3 


0.43 


0.05 


0.001 


CA04 


18.9 


-11.2 


0.72 


0.08 


0.002 


22.9 


-21.3 


2.35 


0.25 


0.009 




25.5 


-15.1 


1.12 


0.15 


0.003 





16.5 


-7.1 


0.39 


0.04 


0.001 


CA05 


18.9 


-8.5 


0.69 


0.08 


0.002 


22.9 


-13.7 


1.97 


0.23 


0.007 




25.5 


-9.4 


1.36 


0.13 


0.004 





16.5 


-4.6 


0.23 


0.04 


0.001 


CA06 


18.9 


-3.7 


0.45 


0.07 


0.001 


22.9 


-3.8 


1.37 


0.19 


0.004 




25.5 


-2.7 


0.75 


0.11 


0.001 



Table 3: The spillover temperature Tg, spillover error Tg >err , zenith opacity r and opacity error T err for 
each antenna and filter. The errors are determined from the fits to derive T$ and r. 



ally last for ~ 15 minutes and show shorter period 
variations (1-10 min), compared to the tens of minute 
timescales associated with clear sky conditions. An- 
tennae 2 and 3 were stowed during this time period, 
and so measure a signal from the zenith, while the 
other 4 antennae were tracking astronomical sources. 
The wavelet transforms for these two antennae are sim- 
ilar and clearly differentiated from the other four, as 
would be expected if they are pointed in different di- 
rections. Furthermore, close inspection of the times 
that each system passes across each antenna shows 
a few minute differences, running from antenna 6 to 
antenna 1 (a distance of 6 km), consistent with wind 
speeds of 10-20 m/s and weather systems arriving from 
the west. Importantly, notably absent from the plots 
is any modulated power which might arise from sys- 
tematic effects, such as thermally induced oscillations 
in the receiver cabin due to the air conditioners. Also 
absent is any signal attributable to RFI. 

In Figure [18] is shown the same analysis, but this 
time for the sensor which measures the cabin tempera- 
ture. Again these plots show similar characteristics to 
one another, but are clearly different to the radiometer 
data. The thermal variations due to the cycling of the 
air conditioners are now clear (and entirely absent in 
the radiometer signal), as are substantial differences 
between the behaviour of the air conditioners in indi- 



vidual antennae. Comparison with Fig. [T7] confirms 
that the measures undertaken to thermally stabilise 
the WVR units have been successful. 

3 Atmospheric Phase Deter- 
mination 

The line strength of the water vapour emission is de- 
scribed in units of brightness temperature Tj, (we also 
use the term "sky temperature" interchangeably) in- 
stead of using specific intensity I v . Tt is defined as 
follows: 

Tb = i 1 ' < 14 > 

where v is the frequency. In the Rayleigh- Jeans 
regime (hv <C kT), the brightness temperature of a 
black body is equal to its physical temperature. 

Below 30 GHz, absorption and emission is dominate d 
by the weak 6i6-5 23 transition of H?Q ijLi ebd ll 19851) . 
iPardo et all (I20011H By follow ing ICarilli et all (| 19991 ) 
as well as llndermuehld l|201lf ) we arrive at the path 
difference (or path excess) Lvx between two lines of 
sight: 
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Figure 13: During calibration in March 2011, sky 
data was recorded and the corresponding sky tem- 
peratures are plotted. The three lines are for the 
PWV values of 5mm, 50mm (the lowest PWV 
value recorded at Narrabri and one order of mag- 
nitude higher) and 23mm (the actual value at the 
time of observation). The profiles were calcu- 
lated using the ATM code. Each filter location 
is marked by three boxes, to show its centre and 
width. 



L VX = 1.763 x 10 -6 APWpw [mm] (15) 

where APWV is the difference in the water vapour 
columns between the two beams in mm, pw is the 
density of water in kg/m 3 and T the temperature in 
K. We thus obtain a path length excess in mm, and 
from this the phase angle difference A$ in degrees, 
with A in mm: 

A$ = 360^^ [deg] (16) 
A 

From equation ll5l we can also establish an approx- 
imate relationship under the assumption of T = 292 
K: 

£vx ~ 6 • APWV [mm] (17) 

We now discuss how to determine APWV from the 
WVRs and so be able to calculate the phase difference. 

4 Extracting Phase 

Phase extraction is achieved by measuring the small 
variations in the line temperatures in each of the four 



Figure 14: Skydip functions for all units. The solid 
line shows the linear fit through the skydip data in 
each filter. The intersect between the fit and the 
a;-axis determines the spillover temperature. The 
highest temperatures are measured at 22.9 GHz, in 
the middle of the water vapour line. Data taken 
on March 7 2011. 

filters at 16.5, 18.9, 22.9 and 25.5 GHz between each 
pair of antennae. Each channel in the WVR is sam- 
pled by the dataset sequentially. This includes all 16 
signal channels including temperature and supply volt- 
age monitoring points, in addition to the tunnel diode 
voltages for the filters. The sampling process currently 
is implemented as follows: 

• For all four signal channels 

— Initial delay of 15 ms to let the input and 
electronic components settle. 

— 8 samples are taken every 5 ms and aver- 
aged over a total timespan of 1,175 ms for 
a total of 1,880 samples. 

• The remaining WVR monitoring points are sam- 
pled with 8 samples only. 

One channel thus takes a total of 1,175 ms to sam- 
ple in this procedure, and the remainder of the mon- 
itoring channels are read at a rate to add to slightly 
less than 5,000 ms total time. Referring to the sensi- 
tivity resulting from the radiometer equation as shown 
in Figure [4] this sampling rate yields a sensitivity of 
about 4 mK for a T rec of ~ 400 K. Refer to Table □ 
for the typical T rec that have been measured. 

4.1 DC Offset 

The measured signal includes unwanted contributions, 
such as the spillover and receiver terms. They are dif- 
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ficult to estimate and moreover may vary with time. 
We work with the differences in signal and assume 
changes in such quantities are slow compared to the 
water vapour induced fluctuations so that they cancel. 
We see in the next section that this indeed allows us 
to follow the phase variations. After converting the 
raw voltages into sky temperatures using the gain cal- 
ibration factors derived from the hot and cold load 
measurements, the total temperature Tt measured in 
each filter can be defined as: 



Tt = Td + Ty + Ts 



(18) 



where Td is the dry air component, Ty is the H2O 
line temperature and Ts is the spillover temperature 
and other instrumental terms. Td will be the same for 
each antenna looking at the same azimuth and eleva- 
tion, whereas Ts, the spillover temperature, is antenna 
specific. The signals for two antennae al and a2 are 
thus 



Tr.ol = Td.ciI + Tv.al + Ts,al 
Tt.u2 = Td,ci2 + Ty,a2 + Ts,a2 



and their difference becomes: 



With To.al = To,a2 



AT T = T V ,al - T V ,a2 + T S ,al ~ T S ,a2 (19) 

Our assumption is that the spillover terms only 
vary slowly. TV will change rapidly on the other hand. 
We can remove slow drifts and offsets by calculating 
differences in a running mean over the time period t 
around each measurement with 2N points: 



T DC {t) = 



t+NAt 

E T D +Tv + Ts 

t-JVAi 

27V 



(20) 



where At is the time step between each sample, 
so that the variation between antennae al and <z2 be- 
comes defined solely by TV: 



ATt — Tv.al — Ty,a2 



(21) 



where these terms are averaged between t — NAt 
and t + NAt. 

4.2 Weighting Coefficients 

From these temperature differences, the water vapour 
excess between two lines of sight can be determined. 
The conversion from temperature differences in each 
filter and antenna pair to electrical path length £y 
is achieved by deriving a weighting coefficient Cw for 
each filter. By modelling the temperatures in the filter 
frequencies and bandwidths for an assumed standard 
atmosphere, we determined the ratio of the filter tem- 
perature and total wet path for each filter to obtain a 
water vapour calibration factor Kf for each filter /: 



Kt 



II 



(22) 



Here £y is the total wet path for the model at- 
mosphere as given by equation 1151 The sensitivity 
to PWV variations is listed in Table 0] and Tf is 



the temperature measured over the filter bandwidth 
for the wet path length, obtai ned by integrating the 
Van V leck - Weisskopf profile (|Tahmoush and Rogers! 
(2000)) over the respective filter frequency and band- 
width. From inspection of Table U it is obvious that 
PWV variations dominate the change in total wet path, 
pressure variations have a small and temperature vari- 
ations a negligible effect. The range of values deter- 
mined for the water vapour calibration factor Kf for 
a variety of atmospheric conditions are listed in Table 
[5] This shows that the calibration factors reflect the 
same sensitivity to PWV already discussed for the wet 
path, however, the relative weights for each filter show 
that the position on the water line affects their weight 
in a given channel. From the calibration factors and 
the relative positions on the water line it is then possi- 
ble to determine the relative weighting coefficients Cw 
for each filter /: 



C W j - 



K7 



K 2 

11 16.! 



l 18.9 

with the normalisation 



7, Cw.f - 1 



(23) 



(24) 



f=i 



The coefficients Cw have been determined for a 
range of atmospheric conditions as listed in Table [6] 
The coefficient with the least weight is at 16.5 GHz, 
which also is the filter least sensitive to atmospheric 
changes, with a weighting less than 2%. The highest 
weight coefficient contributing 60% of the weight is in 
the 22.9 GHz filter band (i.e. near the line centre) and 
varies by ~ 3% from the mean value. 

For the extraction algorithm, the coefficients from Ta- 
ble[S]as determined for an atmosphere of P=1013 hPa, 
PWV=20 mm and T=292 K are used. With calibra- 
tion factors Kf, they are applied to the temperature 
differences obtained from equation [2T] This yields the 
excess path L x j as determined from each filter: 



AT S 
K f 



(25) 



The weighted sum to obtain Cvx is then built with 
/ = filters 1 through 4 by using the weighting coeffi- 
cients derived from the model atmosphere in equation 



/=4 

L V x = ^2 C W J ■ L x>f [mm] (26) 
/=i 

From this it is trivial to obtain the phase differ- 
ence A$ incurred by the path excess Lvx at observing 
wavelength A using equation 1161 

5 Demonstration of Phase Cor- 
rection 

We demonstrate that the WVRs can track phase vari- 
ations by comparing the phase measured directly when 
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Atmosphere 


Variable 


Total wet path L v 


P=1013hPa 


PWV 30 mm 


181.2 mm 


PWV 1 mm 


6.0 mm 


T=292 K 


Spread 


87.6 mm 




1030 hPa 


120.8 mm 


T=292 K 


960 hPa 


120.8 mm 


PWV=20 mm 


Spread 


0.0 mm 


P=1013hPa 


270 K 


130.6 mm 


310 K 


113.9 mm 


PWV=20 mm 


Spread 


8.3 mm 



Table 4: The total wet path L v for each filter. In the first row, PWV is varied while keeping pressure P 
and temperature T constant. The effect ("spread" about the mean) on the total wet path L v is 87.6 mm. 
In the second row P is varied with no effect on L v and in the third row, T is varied, affecting L v by less 
than 10% as compared to varying PWV. 



Water Vapour Calibration Factor Kf [K/mm] 



Atmosphere 


Variable 


16.5 GHz 


18.9 GHz 


22.9 GHz 


25.5 GHz 


P=1013hPa 


PWV 30 mm 


0.044 


0.091 


0.232 


0.164 


PWV 1 mm 


0.043 


0.091 


0.259 


0.169 


T=292K 


Spread 


0.001 


0.000 


0.010 


0.003 


T=292K 


1030 hPa 


0.044 


0.092 


0.239 


0.166 


960 hPa 


0.042 


0.089 


0.248 


0.162 


PWV=20 mm 


Spread 


0.001 


0.001 


0.005 


0.001 


P=1013hPa 


270 K 


0.041 


0.083 


0.207 


0.150 


310 K 


0.045 


0.098 


0.269 


0.177 


PWV=20 mm 


Spread 


0.007 


0.010 


0.031 


0.014 



Table 5: The water vapour calibration factor Kf for each filter and a variety of atmospheric conditions. 
The spread is half the range in path length for the variation in the relevant variable (column 2). The 
largest spread occurs in the 22.9 GHz filter, confirming that this is the most sensitive channel. 









Weighting 


Coefficient Cw 




Atmosphere 


Variable 


16.5 GHz 


18.9 GH 


z 22.9 GHz 


25.5 GHz 


P=1013hPa 


PWV 30 mm 


0.0213 


0.0911 


0.5929 


0.2946 


PWV 1 mm 


0.0172 


0.0791 


0.6353 


0.2685 


T=292K 


Spread 


0.0021 


0.006 


0.0212 


0.0131 


T=292K 


1030 hPa 


0.0205 


0.0886 


0.5998 


0.2912 


960 hPa 


0.0179 


0.0818 


0.6304 


0.2698 


PWV=20 mm 


Spread 


0.0013 


0.0034 


0.0153 


0.0107 


P=1013hPa 


270 K 


0.0228 


0.0930 


0.5798 


0.3045 


310 K 


0.0179 


0.0825 


0.6272 


0.2723 


PWV=20 mm 


Spread 


0.0025 


0.0053 


0.0237 


0.0161 



Table 6: The weighting coefficients Cw for each filter and a variety of atmospheric conditions. The spread 
is half the range in path length for the variation in the relevant variable (column 2). The largest spread 
occurs in the 22.9 GHz filter and amounts to about 3.5%. 
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observing a bright phase calibrator to the phase deter- 
mined from the WVR signals over the same time pe- 
riod. A phase plot from the unresolved astronomical 
calibrator source 0537-441 is compared to that derived 
from water vapour radiometers in Figure QI|] for longest 
baseline of 4500 m formed between antennae 1 and 6. 
We implicitly assume here that the self-calibration un- 
dertaken when measuring the phase calibrator tracks 
the phase perfectly, so that the performance of the 
WVRs can be assessed by comparing the residuals be- 
tween it and the calibrator phases to those between the 
calibrator and the "interpolated" phases. These latter 
phases represent the best estimates that can be ob- 
tained when applying a calibrator phase measurement 
at the beginning and end of an observation sequence, 
and then linearly interpolating between them. In the 
Figure the calibrator phase fluctuations are shown in 
black and the WVR derived fluctuations in red. It 
is evident that they broadly track each other on long 
time scales. On short timescales however, there is con- 
siderable scatter about the true phase values. 

In Figure(2S]the same time span for a short baseline 
of 92 m between antennae 1 and 2 is shown. In this 
case, however, it is not clear that an improvement over 
the interpolated phase has been obtained. However, 
the overall phase noise is also much lower, as expected 
for a short baseline. 

To assess the merit of the WVR corrections we 
compare the standard deviations of the residual phases 
between the WVR and calibrator phase and to the 
residuals between interpolated and calibrator phases. 
For the 4,500 m baseline the residuals are considerably 
less, with a = 14.2° using the WVRs, compared to 
37.2° for the interpolated residual phases. On 92 m 
baseline, however, the residuals have not improved; 
11.5° for the WVRs compared to 9.9° for the inter- 
polated residual phase. 

A second performance indicator we examine com- 
pares the interpolated residual phases with the WVR 
residual phases as a function of baseline length, as 
shown in Figure |2"T1 This shows on long baselines the 
residual phase error being improved by ~ 20° , however 
on short baselines it is worse by ~ 3°. For practical 
purposes this latter degradation is negligible because 
the correlation efficiency e (see equation Sj2]) that cor- 
responds to a 3° higher phase noise is less than 5%. 
This can also be seen from Figure [2] The data do sug- 
gest, however, that a noise floor limitation of about 
10° in phase variation at 22 GHz exists which cannot 
be improved upon using the WVRs. Table [7] lists the 
resulting standard deviations of the residuals and the 
correlation efficiencies for all baselines. As discussed 
above, the correlation efficiency of the WVR residual 
phase (cwvr) is slightly lower than the correlation effi- 
ciency of the interpolated residual phase (ei n t) for the 
short baselines, being reduced by on average 0.02 (but 
remaining in excess of 0.9). On the long baselines, 
however, there is substantial improvement: the cor- 
relation efficiency is improved to above 0.9 when the 
WVR corrections are applied. For the interpolated 
residual phase, the correlation efficiency is only about 
0.63 for the long baselines. 



6 Summary 

A set of water vapour radiometers (WVRs) have been 
developed for the Australia Telescope Compact Array 
that enable the atmospheric phase to be tracked at 
millimetre-wavelengths while the telescope is observ- 
ing a source. These measure the fluctuations of the 
temperature differences in the telescope beam between 
each pair of antennae, through each of four 1 GHz wide 
filters spread across the 22 GHz H2O water vapour line, 
during the 10 s cycle time of the telescope. A set of cal- 
ibration coefficients have been determined between the 
path length difference and the temperature difference 
for each filter based on the properties of a model atmo- 
sphere characteristic of the site conditions at Narrabri. 
A weighted sum of the four temperature differences 
then yields the path length difference between each 
antenna pair, and hence the phase difference at the 
observing frequency. 

The system has been verified by measuring the 
phase variations between antennae pairs recorded on 
a bright phase calibrator and compared to the phases 
determined using the WVRs. When the phase fluctu- 
ations exceed a noise floor (of about 10° at 22 GHz, or 
a path length of ~ 0.4 mm) then the radiometers im- 
prove the phase tracking of the telescope. In practice, 
depending on the weather conditions, this might yield 
improvements on long baselines (where phase varia- 
tions are largest), and slight decorrelations on short 
baselines. It will be up to the user to determine when 
to apply the WVR phase corrections to their astro- 
nomical data. 

Further comparison of the WVR-determined phases 
with phase calibrators under a variety of observing 
conditions is needed to quantify under which condi- 
tions, as a function of baseline and observing frequency, 
the WVRs may be expected to improve the tracking 
of phase during the observation of a source. Such a 
monitoring program may also be able to yield a set 
of empirically derived calibration coefficients for each 
filter, for a variety of weather conditions, that relate 
the path length difference to the measured tempera- 
ture fluctuation. It would be of interest to compare 
these to the atmospheric model-derived coefficients we 
have presented in this paper. 
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Interpolated 


WVR Corrected 


1 V.I V \ V /~\v Tf-) 1 | | , , j | +" i vi 

improvement m 


Baseline [m] 


u Int 


tint 


U WVR 


£\VVR 


correlation efficiency Ae 


92 


9.9 


0.97 


11.5 


0.96 


-0.01 


230 


9.3 


0.97 


11.3 


0.96 


-0.01 


138 
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Spillover temperatures for all antennae and filters 




Figure 15: The spillover temperatures for all an- 
tennas and frequencies. Refer to Table [3] for the 
values. Each filter is shown as a box of 1 GHz 
width and Ts. err height. Additionally, the centre 
point is depicted with a symbol to allow identifi- 
cation of which antenna/ WVR the point belongs 
to. 



Zenith Opacity t far all antennae and filters with ATM Overlays far 
PWV=39, 34,28mm T=300K, P = 969hPa 




16 18 20 22 24 26 



Figure 16: Zenith opacity r as determined from 
the sky dips. Refer to Table [3] for the values. Each 
filter is shown as a box of 1 GHz width and Ten- 
height. Additionally, the centre point is depicted 
with a symbol to allow identification of which an- 
tenna/WVR the point belongs to. Over-plotted 
are three opacities derived from the ATM model 
atmosphere for PWV values of 28mm (lowest opac- 
ity), 34mm (middle line) and 39mm (highest opac- 
ity). 




Figure 17: The Morlet wavelet decomposition of 
a 12 hr time series of raw voltages measured in 
the 25.5 GHz channels on each antenna on 23 May 
2012. Antennae 2 and 3 were stowed, antennae 
1, 4, 5 and 6 were tracking astronomical sources 
during this period. The x-axis is the time of day 
and the y-axis the period (in minutes) for the sig- 
nal power. The solid line marks the cone of influ- 
ence from the edge effects, resulting from the finite 
size of the dataset. Above these lines, the results 
may be influenced by this numerical, though the 
effects are clearly small here. The similarities of 
the transforms between the antennae pointing in 
the same directions are striking, as is the absence 
of any RFI. 
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Figure 20: Phase comparison for data obtained at 
22 GHz on 4 June 2011 on a 92 m long baseline be- 
tween antennae 1 and 2. Calibrator phase - solid 
(black), WVR derived phase - dash-dotted (red). 
The dashed line represents the interpolated phase 
between the measurements and the start the the 
end of the observing period. Calibrator phase mi- 
nus WVR phase - dash-dot-dot (blue; i.e. "WVR 
residual phase"), calibrator phase minus interpo- 
lated phase - dotted (magenta; i.e. "interpolated 
residual phase"). 



Figure 18: The Morlet wavelet transform of the 
same time span as Figure [17] but for the outer shell 
temperature of the WVR units. The activity and 
oscillatory behaviour of the air conditioning units 
in all antennae are clearly evident, as are the dif- 
ferences in the thermal behaviour between receiver 
cabins. 




Figure 19: Phase comparison for data obtained at 
22 GHz on 4 June 2011 on a 4,500 m long base- 
line between antennae 1 and 6. Calibrator phase 
- solid (black), WVR derived phase - dash-dotted 
(red). The dashed line represents the interpolated 
phase between measurements made at the start 
and the the end of the observing period. Calibra- 
tor phase minus WVR phase - dash-dot-dot (blue; 
i.e. "WVR residual phase"), calibrator phase mi- 
nus interpolated phase - dotted (magenta; i.e. "in- 
terpolated residual phase" ) . 







Figure 21: The standard deviation of the WVR 
(A's) and interpolated (+'s) residual phases for 
the data in Figs. [19] and [20] plotted against base- 
line length. Upper plot is for short baselines (up 
to 250 m), lower plot is for long baselines that use 
antenna 6 (i.e. ~ 4 km). 



